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ABSTRACT 

We identify through their X-ray spectra one certain (W37) and two probable (W17 and X4) quiescent 
low-mass X-ray binaries (qLMXBs) containing neutron stars in a long Chandra X-ray exposure of the 
globular cluster 47 Tucanae, in addition to the two previously known qLMXBs. W37's spectrum is 
dominated by a blackbody-like component consistent with radiation from the hydrogen atmosphere 
of a 10 km neutron star. W37's lightcurve shows strong X-ray variability which we attribute to 
variations in its absorbing column depth, and eclipses with a probable 3.087 hour period. For most of 
our exposures, W37's blackbody-like emission (assumed to be from the neutron star surface) is almost 
completely obscured, yet some soft X-rays (of uncertain origin) remain. Two additional candidates, 
W17 and X4, present X-ray spectra dominated by a harder component, fit by a power-law of photon 
index ^1.6-3. An additional soft component is required for both W17 and X4, which can be fit with 
a 10 km hydrogen- atmosphere neutron star model. X4 shows significant variability, which may arise 
from either its power-law or hydrogen-atmosphere spectral component. Both W17 and X4 show rather 
low X-ray luminosities, Lx (0.5-10 keV) ~ 5 x 10 31 ergs s _1 . All three candidate qLMXBs would be 
difficult to identify in other globular clusters, suggesting an additional reservoir of fainter qLMXBs in 
globular clusters that may be of similar numbers as the group of previously identified objects. The 
number of millisecond pulsars inferred to exist in 47 Tuc is less than 10 times larger than the number 
of qLMXBs in 47 Tuc, indicating that for typical inferred lifetimes of 10 and 1 Gyr respectively, their 
birthrates are comparable. 

Subject headings: accretion disks — binaries: close, eclipsing — binaries : X-rays — globular clusters: 
individual (NGC 104) — stars: neutron 



1. INTRODUCTION 

Several soft X-ray transients, identified in outburst 
to be accreting neutron star (NS) systems, have 
been observed in qui escence (see iCamnana et alJ 119981: 
iRutledee et al.l2002b|) . These systems show soft spectra, 
generally consisting of a thermal, blackbody-like com- 
ponent, and in most cases a harder component extend- 
ing to higher energies, usually fit with a power-law of 
photon index 1-2. The nature o f their X-ray emission 
remains unsolved. iBrown et alJ l)1998l see also Cam- 
pana et al. 1998) advanced the idea that the soft ther- 
mal component seen in these systems can be explained 
by the release, over long timescales, of heat injected 
into the deep crust by pycnonuclear reactions driven 
during accretion (the "deep crustal heating" model). 
This scenario accurately predicts the quiescent luminos- 
ity of some qLMXBs. based on their outburst history 
ijRutledge et al]l2000t l2001aD , and must operate at some 
level. (Recent observ ations of low quiescent luminosi- 
ties in some qLMXBs llColpi et al.ll200lUCarnpana et alJ 
12001 lYakovlev et alJl2003t lWiinandsll2004ft may require 
nonstandard cooling processes.) 

However, the deep crustal heating model cannot ex- 
plain the hard power-law component, which is attributed 
to c ontinued accretion and/or a shock from a pulsar win d 
(see iCampana et al.lli"99& iMenou fc McClintockll200l|) . 
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Continued accretion has also been suggested as an ex- 
planation for the thermal component, as the radiation 
spectrum from matter accreting radially onto a neu- 
tron star should be similar to that from the ionized hy- 
drogen atmosphere of a hot non-accreting neutron star 
ijZampieri et al.lll995|) . Models without any continuing 
accretion have difficulty explaining the short-timescale 
(~ 1 4 s) variability observ e d from Aquila X-l an d Cen 
X-4 (|Rut,ledge et al.ll200 2b: Ca mpana et "all 12004?). and 
the days-timescale va riability observed from Cen X-4 
ijCampana et al.l Il997|) and other qLMXBs. However, 
many qLMXBs in globular clusters show neither sig- 
nificant vari ability nor an additi onal power-law spectral 
component ijHeinke et al.ll2003al hereafter HGL03). Un- 
derstanding the emission process will be critical for mod- 
eling the observed spectrum and deriving fundamental 
physi cal parameters, such as mass, radius, and mag netic 
field ijLattimer fc Prakasbll200l IBrown et allll998(l . 

Globular clusters are overabundant in accreting NS 
systems compared to the field, with thirty-e ight probable 
qLMX Bs and active LMXBs known so far l)Heinke et alJ 
2003b). The well-studied globular cluster 47 Tucanae 
(NGC 104; hereafter 47 Tuc) is ideal for X-ray studies of 
its binary populations llGrindla y et al.|| 2001at) due t o its 
close distance (4.85±0.18 kpc. |Gratton et al.ll2003fl and 
low r eddening (E(B-V) = 0.024 ±0.004. lOratton et alJ 
I2003D . Two qLMXBs in 47 Tuc, X5 and X7 (ori ginally 
detected by ROSAT, iVerbunt fc Hasingcr 1998), have 
been spectrally identified by HGL03. HGL03 constrained 
the range of mass and radius space for X7 using several 
different possible assumptions about its spectrum, and 
showed that X5 is an eclipsing system with an 8.67 hour 
period and strong dipping activity. 
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We have shown (HGL03, iHeinke et all l2003bfl that 
the identified qLMXBs in 47 Tuc, and in globular clus- 
ters generally, have little to no hard power-law com- 
ponent in contrast to the best-studied field systems , 
Cen X-4 and Aql X-l l|Rutledge et al.ll2001al |2002b). 
The exceptions thus far have been the transient LMXB 
in NGC 6440, which entered o utburst 13 months af - 
tcr the Chandra observation llin't Zand et alj [2001). 
and t he transient LMXB in Terzan 5 ijWiinandsetaLl 
2003), both relatively bright qLMXBs. If the strength 
of the power-law component (and short- term varia bility) 
measures continuing low-lev el accretion ijCamnana et alJ 
119981 iRutledge et alJl2002a|) . then the absence of these 
indicators in globular cluster systems might be taken 
to indica te extremely l ow lev els of accretion activity. 
Recently, Uonker et alJ l|2004a|) noted that the relative 
strength of the spectral component fit with power-laws 
in quiescent NS transients seems to increase for X-ray lu- 
minosities significantly larger or smaller than Lx ~ 10 33 
ergs s _1 , wher e its relative str ength is often at a min- 
imum (but cf. iWiinands et al.l [20031. They suggested 
that this spectral component arises from different pro- 
cesses (accretion vs. pulsar wind shock?) in systems with 
X-ray luminosities above or below 10 33 ergs s _1 . 

In late 2 002, we obtained de eper Chandra observations 
of 47 Tuc (Hei nke et al.ll2004l) . allowing detailed spectral 
modeling of X-ray sources as faint as Lx ~ 10 31 ergs s . 
Most sources were found to be consistent with powerlaw 
or thermal plasma models, but some sources required 
more complicated models. The highly variable X-ray 
source W37 (=CXOGLB J002404.9-720451) showed a 
steep spectrum well-described by a blackbody, and the X- 
ray sources W17 (=CXOGLB J002408. 3-720431) and X4 
(=W125, or CXOGLB J002353.9-720350) showed two- 
component spectra in which the soft component could 
be well-represented by a model of a hydrogen atmo- 
sphere appropriate for a neu tron star. X4 was first iden - 
tified in ROSAT HRI data l)Verbunt fc Hasingerlll99^ . 
while W17 and W37 are located in areas too crowded to 
have been identified in lower-resolution ROSAT data at 
their current fluxes. None of these three have been posi- 
tively i dentified with opti cal counterparts of any category 
l)Edmonds et al.ll2003alH ). In this paper, we study these 
three strong qLMXB candidates in detail. We describe 
the observations and our reduction in § EI the timing 
analysis in § 12.11 and the spectral analyses in § 12.21 Dis- 
cussion and conclusions are in § 

2. OBSERVATIONS AND ANALYSIS 

The data used in this paper are from the 2000 and 
2002 Chandra observations of the globular cluster 47 Tuc. 
Both sets of observati ons and their initial reduction are 
described in detail in IHeinke et alJ (120041): prior analy- 
ses of the 20 00 dataset are described in iGrindlav et aD 
(|2001aU2002h and HGL03. Both observations interleaved 
short exposures using subarrays (chosen to reduce pileup, 
which is not a significant issue for the three sources in 
this paper) with longer exposures. The observations are 
summarized in Tabled we note that OBSJD 3385 suf- 
fered strongly increased background, so we do not use 
that short observation. We reprocessed (using CIAO 3.0) 
both the 2000 and 2002 observations to remove the C/25 
pixel randomization added in standard processing, im- 
plement the CTI correction algorithm on the 2000 data, 



and use updated (Feb. 2004, CALDB v. 2.26) calibra- 
tion files. We did not remove events flagged as cosmic-ray 
afterglows by the standard processing. 290 sources were 
detected using the WAVDETECT routine in the 0.3-6 
and 0.3-2 keV energy bands, and ten additional sources 
(missed by WAVDETECT due to crowding) were added 
to the source list by hand. 

We used the ACIS_EXTRACT routine (|Broos et alJ 
2002) from the Tools for X-ray Analysis (TARA 4 ) to 
construct polygonal extraction regions generally chosen 
to match the 90% encircled energy (at 1.5 keV) from 
the Chandra point-spread function (PSF). We extracted 
source spectra and background spectra from nearby 
source-free regions, and corrected the ancillary response 
function (ARF) for the energy-varying fraction of the 
PSF enclosed by the extraction region. We used the 
Jan. 2004 release of the ACIS contamination model to 
correct for the hydro carbon buildup on the detectors 
({Marshall et alJl2004j) . We adjusted all event times to 
the solar system barycenter using satellite orbit files pro- 
vided by the Chandra X-ray Center. 

2.1. Timing Analysis 

Kolmogorov-Smirnov (KS) tests on the 2002 event ar- 
rival times from W37 and X4 confirm that both are 
clearly variable on short timescales, with KS probabil- 
ities of constancy reaching 10 -47 and 3 x 10 -4 for the 
first 2002 observations of W37 and X4 respectively. W17 
shows only marginal evidence for variability, with KS 
probabilities reaching only 0.04. X4 falls on a chip gap 
in the 2000 data, and W37 is much fainter in the 2000 
data. W37 receives a photon flux over five times higher in 
the first 2002 observation than in any other observation, 
and 75% of the 2002 photons from W37 arrive in the first 
quarter of the data. X4 appears to have dimmed slightly 
between the first and fourth long 2002 observations (at 
99% conf.). 

We show the first (and most variable) sections of the 
2002 lightcurves for W37 and X4 in Figures □ and 
Generating a power spectrum for W37 (using XRONOS, 
5 ) gives a possible peak around 11000 seconds along with 
strong red noise. An epoch folding search using all the 
2002 observations finds a periodicity at 11112.5 (±15) 
seconds. Careful inspection of W37's lightcurves reveals 
likely eclipses, indicated in Figuresn]and|31(left; a folded 
lightcurve), consistent with the 11112.5 second period. 
The two clearest eclipses are separated by 6.2 hours, so 
we also include a lightcurve folded on twice our pre- 
ferred period (Figure |3| right), showing two apparent 
eclipses. This indicates that our preferred period is prob- 
ably correct, but the low count rate of the data may allow 
for alternative solutions. The folded lightcurve suggests 
eclipses 550 ± 200 seconds in length. 

We note that W37's lightcurve is similar to that of 
X5, presented in HGL03 and Heinke et al. (2005, in 
prep), suggesting that the strong variability seen in W37 
is also probably due to a varying column of cold gas 
associated with an accretion disk viewed nearly edge- 
on. However, W37 may show residual emission dur- 
ing eclipses, which is not seen for X5 (HGL03). The 
period indicates a companion mass of order 0.34 Mq, 

4 http:/ /www. astro. psu.edu/xray/docs/TARA/ 

5 http:/ /heasarc. gsfc.nasa.gov/docs/xanadu/xronos/xronos. html 
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if the companion ha s a density appro priate for lower 
main sequence stars ijFrank et al.lll992l eg. 4.9). Since 
these companions are sometimes bloated (|KinjL_fe_ilittei 
| 1999f [Podsiadlowski et all Eool lOrosz fe van Kerkwuk 
20031 iKaluznv fc Thompson! 12003). the true mass may 
be smaller. Assuming a 1.4 M Q neutron star, the orbital 
separation would be ~ 9 x 10 10 cm. 

We calculate the hardness ratio for W37 as the ratio of 
the counts detected in the 1-6 keV and 0.3-1 keV bands 
(this choice divides the detected counts roughly evenly), 
and plot this ratio in the lower panel of Figure ^ Com- 
parison with the upper panel shows that decreases in the 
observed flux correlate well with increasing hardness. We 
confirm this with a Spearman rank-order correlation test 
ijPress et al.lll992 j) comparing the binned (1200 second) 
lightcurve with the similarly binned hardness ratio. This 
gives r s =-0.446, a negative correlation, with chance prob- 
ability 5.6 x 10~ 4 . This is a strong sign that the vari- 
ation is due to changing photoelectric absorption. We 
divide W37's long 2002 observations into high and low 
flux states (e.g., the high flux state portions indicated in 
Figure nj, which we use to extract spectra. 

X4's total and hardness ratio lightcurve for OBSJDs 
2735 and 3384 are shown in Figure [3 No periodicities 
are apparent in X4's lightcurve or power spectra. X4's 
spectrum appears to harden as it brightens, in contrast 
to W37, which indicates the variability is not due to in- 
creasing photoelectric absorption. We confirm this with 
a Spearman rank-order correlation test on OBSJD 2735, 
finding r s =0.444 and chance probability 6. Ox 10 -4 for the 
binned (1200 second) lightcurve and hardness ratio. We 
extract high and low portions of X4's lightcurve from the 
first observation for spectral fitting, which are marked in 
Figure [3 

2.2. Spectral Analysis 

In lHeinke etafl p004j) . we identify W37, X4, and W17 
as possible qLMXBs by their agreement with models con- 
sisting of a hydrogen-atmosphere model and a power-law, 
along with their failure to fit single-component mod- 
els typical of thermal plasma (such as seen in cata- 
clysmic variables and X-ray active binaries). Here, we 
analyze time-resolved spectra for W37 and X4, and use 
the merged spectra (2000 and 2002) for W17 (since it 
shows no variability). We use the XSPEC phabs model 
to des cribe neutral absorbing gas. We use the Lloyd 
(2003) model to describe thermal blackbody-like emis- 
sion from the hydrogen atmosphere of a neutron star, 
fixing the redshift to 0.306 (appropriate for a 1.4 M©, 10 
km neutron star) an d take the distance to be 4.85 kpc 
ijGratton et al.ll2~003D . The total number of counts in the 
2002 spectra of W37, W17 and X4 were 1277, 909, and 
1426, while from the 2000 data 46, 114, and 62 counts 
were available. We binned spectra with more than 1000 
counts with 40 counts/bin, spectra from the 2000 obser- 
vations or with fewer than 50 counts at 10 counts/bin, 
and other spectra at 20 counts/bin. 

2.2.1. W17 

W17 shows no variability on short or long timescales, 
so we fit the combined 2002 and combined 2000 datasets 
simultaneously. We try a variety of absorbed single- 
component models including blackbody, bremsstrahlung, 
MEKAL, power-law, and hydrogen atmosphere neutron 



star models, but do not find good fits with any. The 
best fit of these, a power-law model, gives \l = 1-36, for 
a 2% null hypothesis probability (nhp). Models consist- 
ing of two mekal components also do not produce good 
fits. Adding a soft component, s uch as the hy drogen- 
atmosphere neutron star model of iLlovdl l|2003f) . to the 
power-law model improves the statistics tremendously 
(xl = 1-05, nhp=37%; F-test gives probability of 9 xl0~ 5 
that the extra component is not needed). The inferred 
radius of the neutron star model, R — 11.8jij'g km, is 
perfectly consistent with the canonical 1.4 M©, 10 km 
neutron star predictions, although it cannot usefully con- 
strain the neutron star structure. (The inferred radius 
for a blackbody would be only 1.2±°-\ km.) A qLMXB 
model is thus an excellent explanation for this spectrum, 
shown in Figure The parameters for a fit with a 
Lloyd H-atmosphere model and power-law are listed in 
Table The power-law component provides a majority 
(65±22%) of the 0.5-10 keV flux, and the thermal com- 
ponent is quite weak, with unabsorbed X-ray luminosity 
Lx(0.5 - 2.5) = 1.7 x 10 31 ergs s _1 and inferred bolo- 
metric luminosity L^ i = 6.7 x 10 31 ergs s _1 . 

2.2.2. W37 

We identify pieces of W37's lightcurves at high and low 
flux levels in the four long observations, with only a low 
flux level in the second long observation; the (barycen- 
tered) times are listed in Table El We extracted spectra 
from these periods, as well as from the summed 2000 
observations, and background spectra from surround- 
ing source-free regions, producing eight spectra of varied 
quality. 

The hardening of the spectrum with decreasing flux 
suggests that variation of the column density is responsi- 
ble for the spectral variations, and so we allow the gas col- 
umn to vary between spectral segments. However, no ab- 
sorbed single-component model was able to describe the 
eight spectra simultaneously by only varying the absorp- 
tion and one other model parameter. (Allowing all the 
parameters to vary freely allows some reasonable fits, but 
loses physical meaning. For instance, a blackbody with 
strongly varying normalization, temperature and absorb- 
ing column, while a reasonable fit, is difficult to inter- 
pret.) No single- or double-temperature thermal plasma 
model (represented by VMEKAL models in XSPEC, set 
to cluster abundances) provides a good fit to the bright- 
est portions of the data. The brightest portions of the 
data are best fit by a blackbody or hydrogen atmosphere 
neutron star model with varying absorption. If described 
by a blackbody, the inferred radius is only 1.4 ± 0.2 km 
in size. This component is well-described by a hydrogen- 
atmosphere neutron star model with a constant physical 
radius of order 10 km (12. 3^3 g km) and constant tem- 
perature (82±* eV). Therefore we infer that W37 is a 
qLMXB containing a neutron star. 

The fainter portions show additional soft flux beyond 
what can be accounted for by an absorbed H-atmosphere 
(or blackbody) component. This suggests a second com- 
ponent that does not suffer the same absorption as the 
primary spectral component. We model this component 
with a separately absorbed power-law, although a vari- 
ety of models can account for this low-count component 
equally well (including a 0.6 keV bremsstrahlung spec- 
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trum, a 0.1 keV blackbody with inferred 1 km radius, 
or a 0.13 keV thermal plasma model), and we do not 
know its physical origin. We hold the faint component, 
its absorption, and the H-atmosphere model parameters 
fixed between data portions, and vary only the absorp- 
tion to the H-atmosphere model. This is the XSPEC 
model PHABS ( POW + PHABS * HATM ), varying 
only the second PHABS component between observa- 
tions. We find a good fit (xl — 1-24 for 57 degrees of 
freedom, or dof, and an nhp of 10.6%). We are unable to 
impose serious constraints upon this model if additional 
parameters are allowed to vary, since the faint portions 
of the data have very few counts. Therefore we cannot 
usefully constrain any temporal or spectral variability in 
the faint component. We show part of this fit (only 4 
of 8 spectra to reduce confusion) in Figure and give 
details of the fit in Tables [21 and |SJ The fitted absorption 
column to the probable neutron star varies by a factor 
> 100, decreasing to a value consistent with the cluster 
column during parts of the first 2002 observation. The 
faint component is responsible for a very small portion 
(2~t\ %) of the unabsorbed flux, and produces a photon 
index of 3j^. The absorption column to the faint com- 
ponent is poorly constrained (Nh = 4.5^4 5 x 10 20 ) but 
consistent with the column to the cluster. The lowest- 
flux spectrum (from OBSJD 2736), dominated by the 
faint component, is one of the spectra shown in Figure^ 

2.2.3. X4 

X4 varies during the first 2002 observation, but is not 
clearly variable within other observations. We separately 
fit the high-flux and low-flux parts of the first 2002 obser- 
vation (OBSJD 2735), and data from each of the other 
three long 2002 observations, plus one spectrum from 
the combined 2000 observations. We found that the re- 
maining 2002 observations were spectrally indistinguish- 
able, so we combined all 2002 data except OBSJD 2735 
into one spectrum, and fit the high and low-flux parts 
of 2735, the rest of the 2002 data, and the combined 
2000 data simultaneously. Fits with any single absorbed 
component, with all parameters allowed to vary, failed 
to give reasonable fits to the spectra. A power-law fit 
gave the lowest xl = 1-45, nhp=3.9%, if Nh is allowed 
to vary substantially along with the power-law parame- 
ters. However, if Nh is not allowed to vary during the 
2002 data (as indicated in § 12.11 above - ), then x 2 =1.52, 
nhp=1.9%. Adding a H-atmosphere model improves the 
fit substantially (F-statistic=4.38, prob. 1.9% extra com- 
ponent not required), and results in an inferred radius 
consistent with the 10 km canonical NS radius. We con- 
clude that X4 is a good candidate for a qLMXB. Alter- 
nate models involving two VMEKAL components with 
all temperatures and normalizations free do not produce 
good fits. 

The variability during OBSJD 2735 demands that at 
least one component of the spectrum varies. Varying 
Nh alone does not produce a good fit (xl = 2.41, 
nhp=l x 10 _4 %). Varying the power-law component docs 
produce a good fit (xl = 1-30, nhp=ll%; see Table©- 
We show this fit for the two components of OBSJD 
2735 (Fig. and for the remainder of the 2002 obser- 
vations, indicating the contributions of the two compo- 
nents (Fig. EI). Varying only the H-atmosphere com- 



ponent's temperature and radius also produces a good 
fit (xl = 1-35, nhp=7.3%; see Table©. Varying both 
the power-law and H-atmosphere parameters (xl = 1.21, 
nhp=20%) is only a marginal improvement over varying 
a single component, and requires correlated swings in nu- 
merous parameters (which we assume is less likely than 
one or two varying parameters). Varying only the power- 
law normalization also produces a good fit (xl = 1.33, 
nhp=7.5%), with a power-law photon index of 2.2^q'3. 
Varying the H-atmosphere temperature alone does not 
produce a good fit (xl = 1-62, nhp=0.7%), nor does 
varying the normalization (and inferred radius) alone 
(xl = 1.85, nhp=0.07%). If the thermal component 
is due to continued accretion, and the variation is due 
to changes in the accretion rate onto the neutron star, 
we would expect only the temperature to vary, not the 
inferred radius. It is possible that the emitting area 
varies if accretion is ongoing, but variation in the emit- 
ting area has not been clearly seen in fie l d qLMXBs with 
highe r statistics (jRutledge et al.ll2002ri iCanroana et alJ 
I2004|) . Therefore we conclude that the power-law prob- 
ably varies, and that we do not have strong evidence for 
variability of the thermal component in this source. 

We can also compare X4's flux in our Chandra ob- 
servations to its flux in ROSAT HPJ observations from 
1992-1 996 (which contain very lit tle spectral informa- 
tion). iVerbunt fc Hasingeil 1)19981) report 49±9 counts 
from X4 (which we confidently identify with our source) 
in 58820 seconds of exposure. Using a power-law spec- 
tral fit to only the 0.5-2.5 keV Chandra data (with Nh 
fixed to the cluster value, 1.3 x 10 20 cm -2 ), we find a 
power-law index of 2.9±0.1, with (absorbed) 0.5-2.5 keV 
X-ray fluxes ranging between 2.7 x 10~ 14 to 1.0 x 10~ 14 
ergs s _1 cm -2 (for different parts of the Chandra data, 
as above). Using this spectrum, the PIMMS tool 6 gives 
an absorbed 0.5-2.5 keV X-ray flux of 2.8 x 10 -14 ergs 
s _1 cm" 2 for the ROSAT HRI data. This flux is consis- 
tent with the upper end of the range of fluxes observed 
with Chandra . 

3. DISCUSSION 

The 47 Tuc X-ray sources W37, W17 and X4 are very 
probably qLMXBs. No other identified source types 
observed in globular clusters show spectral components 
consistent with the emission from a 10 km radius hy- 
drogen atmosphere, as observed in these objects. Mil- 
lisecond pulsars show thermal emission f rom a much 
small er area, and a t lower luminosities llZavlin et alJ 
l2002tlGrindlav et al.ll20?)l iBogdanov et alJl2004D . Mag- 
netic CVs show soft, blackbody- like components plus 
harder emission, but with emission radii that are much 
larger (generally hundreds of km) than the implied radii 
(0.5-2 km) from blackbody fits to these sources. We also 
note that these three objects are the three X-ray bright- 
est objects (excepting X7, another known qLMXB) not 
yet positively identified in our HST o ptical identifica- 
tion program (Ed monds et al.ll2003albf) . If W17 were a 
cataclysmic variable (or, even less likely, an X-ray ac- 
tive binary), its faint upper limit of U > 24 would pro- 
duce an unusua l X-ray to optical flux ratio, which led 
lEdmonds et al.l (|200 3a*) to suggest this object as a possi- 
ble qLMXB. X4 fell outside the HST fields of view ana- 

6 Available at http://asc.harvard.edu/toolkit/pimms.jsp 
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lyzed bv lEdmonds et all l)2003a|) . but its X -ray spectral 
similarities to W17 are convincing. W37's error circle 
showed no blue or variable objects, but two bright main- 
sequence stars fall within a 2a error circle (not unusual 
considering its projected location at the center of the 
cluster) , likely obscuring the t rue (fainter) companion (as 
for X7, Edmonds et alJl2002fl. T he far-UV imaging ob- 
servations of lKnigge et al.l |2002) include W37's position, 
and may be able to provide stronger constraints on a pos- 
sible UV counterpart. We conclude that W37's varying 
Nh and eclipses indicate that it is a transiently accret- 
ing qLMXB, with marked similarities to X5 (HGL03). 
We also note that W37 is the second qLMXB without a 
recorded outburst to have an identified period, after X5, 
and that its short period implies a very dim low-mass 
companion. 

We can not rule out possibilities other than qLMXBs 
for W17 and X4. It is possible that W17 and X4 
are neutron stars without Roche-lobe filling low-mass 
companions, accreting fro m the intracluster medium 
( Pfah l fe Rappaportl l2001[) or from the normal stellar 
winds of low-mass stellar co mpanions in close orbi ts, in 
so called pre-LMXB systems l|Willems fc Kolbl2003fl . Ei- 
ther of these is possible, but we judge a stellar wind 
model to be more likely, because neutron stars that ac- 
crete from the intracluster medium must also be able to 
accrete the higher- density wind from a low-mass star, 
and because all the (slowly spinning) neutron stars in 
the cluster should display signals of accretion from the 
intracluster medium if any do. The low bolometric lumi- 
nosities from the neutron star surfaces of W17 and X4 
indicate a very low level of t ime-averaged mass trans- 
fer in the iBrown et alJ l)1998fl model, or enhanced neu- 
trino cooling. Mass transfer from a pre-LMXB by a 
stellar wind thus seems a reasonable possibility. How- 
ever, similarly low quiescent luminosities and power-law 
dominated spectra have recently been found for sev- 
eral transient LMXBs: L x (0.5-10 keV )=5 x 10 31 ergs 
s" 1 fo r SAX J1808. 4-3658 in quiescence l)Campana et al.l 
1200^ . L^(0.5-10 k eV)=9 x 10 31 ergs s " 1 (d/8.5 kpc) 2 
forXTE J2123-058 (|Tomsick et al.l2004|) . and L x (0.5-10 
keV)=l x 10 32 ergs s~ 1 (d/4.9 kpc) 2 for SAX J1810.8-2609 
(|.Tonker et al.ll2004bj) . This suggests that W17 and X4 
may be normal qLMXBs. They seem to be generally con- 
sistent with the anticorrelation between X-ray luminos- 
ity, and the fraction of that lu minosity seen i n the p ower- 
law component, suggested bv l.Tonker et alJ l)2004al) . On 
the other hand, W37 (like the qLMXB in NGC 6397, 
iGrindlav et al.ll2001b|K appears to show a smaller power- 
law component than expected from the anticorrelation of 
I.Tonker et all l|2004aft . 

The faint soft component in W37's spectrum must be 
associated with W37, since the position of the source 
does not change between the brightest and faintest in- 
tervals (as would be expected if another X-ray source 
produced this component). The X-rays might be pro- 
duced by th e rapidly rotating companio n star's corona, 
as sugge sted Bildst en fc Rutledgel l|2000D and disproved 
l)Lasotall2000t iKong et al.ll2002j) for quiescent black hole 
systems. We can estimate the coronal saturation lumi- 
nosity of a companion as logL cor = — 2.9 + 21ogi?£(, where 
L cor and Rd are measur ed in units of the sola r luminosity 
and radius respectively ( Flemin g et alJfl989|) . Assuming 
a companion mass of 0.3 Mq (appropriate for a lower 



main-sequence star in a 3-hour orbit; see § I2.1[1 . we com- 
pute a radius 2.5 x 10 10 cm, and thus a coronal saturation 
luminosity \ogL cor < 29.8 ergs s , which is significantly 
less than the 5 x 10 30 ergs s _1 required for the addi- 
tional component. Scattering of W37's X-rays in an ac- 
cretion disk wind might produce a similar spectrum, as 
obser ved in the eclipsing CVs OY Car DQ Her and UX 
UMa ijNavlor et al.lll988t IMukai et alJl2003l iPratt et alJ 
l2004[) . We are not able to make predictions about the 
efficiency of scattering in such a wind, but we do not see 
any outstanding objections to this scenario. 

Alternatively, this faint less-absorbed component 
might be identified with the power-law component ob- 
served in many field qLMXB s (e.g., iCampana et alJ 
fl998t iRutledge et al.ll200Tblal) . and in W17 and X4. If 
this identification is correct, then this component must 
be generated in an extended environment around the 
neutron star, and not near the surface. A rough esti- 
mate of its size can be estimated by assuming that the 
place of its generation is large compared to the size of 
variations in the accretion disk. (This assumes that the 
variations in Nh observed for the neutron star compo- 
nent are due to material in the accretion disk.) The 
accretion disk is probably of order 4 x 10 10 cm in ra- 
dius, for typical masses of the components (1.4 and 0.1 
Mq) and an accretion disk size one- half the orbital sep- 
aration. Assuming that variations in the accretion disk 
height may be up to one-tenth the radius of the disk, 
this suggests that the environment that generates the 
faint soft component should be larger than 4 x 10 9 cm. 
This is significantly larger than the light-cylinder radius 
for neutron stars spinning at millisecond periods. If the 
X-rays are generated within or at the edge of the mag- 
netosphere, the size of the emitting region implies that 
a pulsar mechanism is operating, and that the observed 
X-rays may be a shock f rom t he pulsar wind as sug- 
gested bv lCampana et all (^998). The physical nature of 
the additional component to W37's spectrum cannot be 
conclusively determined from this small amount of data, 
but the possibility that we are seeing a pulsar wind in a 
system with an accretion disk is intriguing. 

W17 and X4 indicate that a class of faint qLMXBs 
dominated by harder nonthermal emission may exist 
in similar numbers as the class of qLMXBs domi- 
nated by_^hermal e^nission from the neutron star sur- 
face l|Heinke et alJ l2003b[) . Additional qLMXBs with 
quiescent properties similar to SAX J1808. 4-3658 hid- 
ing among the low-luminosity sources in 47 Tuc can- 
not be ruled out, since the ne utron star atmospher e 
is not always clearly detectable ijCampana et al J 12002). 
These particular identifications would be difficult with 
shorter exposures, as have been performed for other 
globular clusters, due to the burying of their telltale 
thermal components under a stronger nonthermal com- 
ponent and/or a blanket of absorbing gas, and their 
significantly lower fluxes. However, it will be worth- 
while to look for similar objects in the nearest, best- 
studied clusters. The numbers of faint qLMXBs with 
Lx < 10 32 ergs s _1 may be similar to those of the 
identified qLMXBs. We note that there are few hard 
qLMXBs among sources with Lx — 10 32 ~ 33 ergs s^ 1 , 
using the results of detailed X-ray and optical studies 
of 47 Tuc, NGC 6397, and NGC 6752 llEdmonds'et~aTl 
I2003albt loTrndlay et a,l.ll2001l4 IPoolev et al.M2002fl. hut 
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a significant population below Lx < 10 32 ergs s" 1 can- 
not be ruled out. If we extrapolate from these two faint, 
nonthermally-dominated qLMXBs to the globular clus- 
te r system (using thei r close encounter frequencies, as 
in lHeinke et alJl2003rA and Poisson statistics), an addi- 
tional 56^36 such systems are sugg ested, added to the 
~100 qLMXBs alre ady anticipated l|Poolev et all 12003 
iHeinke et aT]l2003b|) . Study of faint qLMXB systems 
may shed light on the interior physics of neutron stars, 
the physics of accretion at low mass transfer rates, and 
possibly the tran sition as LMXBs eyqlve to become mil- 
lisecond pulsars l)Burderi et alJl2002t 1200^ . 

Early estimates of the numbers of millisecond pul- 
sars (^10000) and their suggested progenitors, LMXBs 
(12), in globular clusters indicated (for LMXB life- 
times 1/10 those of millisecond pulsars, ~1 and 10 
Gyr) a large discrepancy in their resp ective birthrates 
UKulkarm et a.U fiwl ITTut et al] ITflfll^. thus implying 
that other formation pr ocesses may have created some 
globular cluster pulsars l)Bailvn & Grindlavl lT990). Ra- 
dio timing surveys of 47 Tuc have identified 22 millisec- 



ond pulsars ijCamilo et alJl2000t IPTeire et al]l2001[K and 

X-ray, optical and radio imaging surveys indicate that 
the total nu mber of millisecond pulsars in 47 Tuc is 
of order 30 iHeinke et all 120041: lEdmonds et aD I2003U: 
IMcConnell et al.l 120041 IGrindlav et al.l 12002ft . When 
these constraints are compared with the current number 
of likely quiescent LMXBs in 47 Tuc (5), the birthrate 
discrepancy between millisecond pulsars and LMXBs dis- 
appears, indicating that LMXBs are probably sufficient 
to produce the observed millisecond pulsars. 
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Orbital Phase 



Fig. 1 — Top: Lightcurve of W37 (0.3-8 keV), in 300 second bins, from first two 2002 observations (OBSJDs 2735 and 3384, 65 and 5 
ksec, respectively). Time is labeled in units of the best-fit period, 11112.5 seconds, with integer phases occurring at the times of mid-cclipsc. 
Horizontal bars indicate portions of the data taken for a high-state spectrum (2735, H); the remainder of OBSJD 2735 is taken as a low-state 
spectrum (2735, L). Bottom: Hardness ratio (1-6 keV/0.3-1 keV) lightcurve for W37, corresponding to the same timespan as above, but 
with 1200 second bins. An anticorrclation can be seen between hardness and count rate. 
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Fig. 2. — Top: 0.3-8 keV lightcurve of X4 (W125), in 1200 second bins, from first two 2002 observations. Time is labeled in hours. A 
horizontal bar marks the portion of the data taken for the high-state spectrum. Bottom: Hardness ratio (1-6 keV/0.3-1 keV) lightcurve for 
X4, with 5000 second bins. A correlation can be seen between hardness and count rate. 




Phase Phase 

Fig. 3. — Left: Folded lightcurve of all W37 2002 data on 11112.5 second period, with the data repeated over two phases for clarity. A 
clear eclipse can be seen at phase 0. Each bin is 185.2 seconds long. Right: Folded lightcurve of all W37 2002 data on 22225 second period, 
showing eclipses at phases and 0.5. 
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0.1 1 10 

channel energy (keV) 

Fig. 4. — X-ray spectrum of W17 from the 2002 observations, fit with an absorbed hydrogen atmosphere model and power-law as in 
Table l2l The contributions from the two components are shown, with the power-law dominating at higher energies. 




0.1 1 10 

channel energy (keV) 

Fig. 5. — X-ray spectra of W37, taken from four different portions of the data with different fluxes. The spectra were fit with an absorbed 
hydrogen atmosphere model plus a separately absorbed power-law. Only the absorption on the hydrogen-atmosphere model was allowed 
to change between data segments. Eight portions were fit to produce the results in Tables [2] and [3j but only four are shown here to reduce 
confusion. From highest to lowest flux, these are [2735, H], dashed line; [2735, L], squares and dotted (red) line; [2737, L], triangles and solid 
(blue) line; [2736], stars and dash-dotted (green) line. (See the electronic edition of the Journal for a color version of this figure.) 
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Fig. 6. — X-ray spectra of X4 from two parts of the first 2002 observation, fit with an absorbed hydrogen atmosphere model and a 
power-law (with the power-law component allowed to vary between observations) as in Table |21 The lower-flux spectrum (red) shows the 
greatest differences at higher energies, indicating the difference is not due to photoelectric absorption. (See the electronic edition of the 
Journal for a color version of this figure.) 




Fig. 7. — X-ray spectra of X4 from the rest of the 2002 observations (except those in Fig. [fjj, fit with an absorbed hydrogen atmosphere 
model and a power-law (with the power-law component allowed to vary between observations) as in Table |5] The contributions from the 
two components are shown, with the power-law dominating at higher energies. 



TABLE 1 

Summary of Chandra Observations 



Seq,OBS.ID 


Start Time 


Exposure 


Aimpoint 


Framctime 


CCDs 


300003,078 


onnn a /r i a n7 10 on 

2000 Mar lb 07:18:30 


3875 


A PIG T 

ACIS-I 


0.94 


1 i A 

1/4 


300028,953 


2000 Mar 16 08:39:44 


31421 


ACIS-I 


3.24 


6 


300029^954 


2000 Mar 16 18:03:03 


845 


ACIS-I 


0.54 


1/8 


300030,955 


2000 Mar 16 18:33:03 


31354 


ACIS-I 


3.24 


6 


300031,956 


2000 Mar 17 03:56:23 


4656 


ACIS-I 


0.94 


1/4 


400215,2735 


2002 Sep 29 16:59:00 


65237 


ACIS-S 


3.14 


5 


400215,3384 


2002 Sep 30 11:38:22 


5307 


ACIS-S 


0.84 


1/4 


400216,2736 


2002 Sep 30 13:25:32 


65243 


ACIS-S 


3.14 


5 


400216,3385 


2002 Oct 01 08:13:32 


5307 


ACIS-S 


0.84 


1/4 


400217,2737 


2002 Oct 02 18:51:10 


65243 


ACIS-S 


3.14 


5 


400217,3386 


2002 Oct 03 13:38:21 


5545 


ACIS-S 


0.84 


1/4 


400218,2738 


2002 Oct 11 01:42:59 


68771 


ACIS-S 


3.14 


5 


400218,3387 


2002 Oct 11 21:23:12 


5735 


ACIS-S 


0.84 


1/4 



Note. — 



Times in seconds. Subarrays are indicated by fractional numbers of CCDs. 



TABLE 2 
X-ray Spectral Model Parameters 



Obs 


F x , abs (0.5-10 keV) 
(10~ 14 ergs s" 1 ) 


N H 
(10 20 cm" 2 ) 


kT 

(eV) 


R Lx,jvs(0.5-10 keV) L boLNS 
(km) (10 31 ergs s" 1 ) (10 31 ergs s" 1 ) 


r 


^x,pl(0.5-10 keV) 
(10 31 ergs s" 1 ) 


X 2 /DoF(nhp) 




W17, 2000 and 2002 data 




i tH-U-3 
1 -°~0.4 


*' x -1.3 




15.0+5 5 3 1.5 7.3 


1 q+U-4 


3.3 


1.28/43(10%) 


W37, 2000 and 2002 data 




a 


a 


82 +i , u 


12.3+3g 26 53 


o n +4.8 
a,u -0.9 


0.56 


1.24/57(11%) 










and H-atmosphere fixed, power-law varies 








2735 high 
2735 low 
Other 2002 


q 7+1.3 

1 q+0.5 
1,o -0.7 

2 0+ ' 4 


5.0^ 


53±^ 


10.8+47 2 -° 7 - 3 


2 4+ U ' 3 
o ,+0.5 

n q+0.4 
z -°-0.3 


9.8 
2.5 
4.6 


1.30/38(11%) 


2000 


i q+0.7 








2 4+ 10 

z -^-1.0 


2.4 




X4, Nh and power-law fixed, H-atmosphere varies 



2735 high 


q 9+U.3 

°- z -2.0 
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1.35/38(7%) 


2735 low 


i 7+0.3 
1.1 






19+ 22 


2.0 
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Other 2002 


i.9±S:? 




mtll 


- n +2.0 
'• u -2.6 


2.7 


7.1 








2000 


1.6+1;? 




28+ 32 b 


' u -62 


1.5 


24 









2 
s. 

CD 
-J. 
r. 

§ 



Note. — Spectral fits to W37, W17 and separate parts of X4 data, using Llovd I 2003) hydrogen-atmosphere neutron star model plus a power-law model. For X4, some 
parameters are held fixed between observations, while others arc allowed to vary. All errors are 90% confidence limits. Distance of 4.85 kpc is assumed. Flux measurements 
are absorbed (no correction for iVfj), while Lx and L^ a i arc unabsorbed. Neutron star H-atmosphcrc radius and temperature for assumed grav. redshift of 0.306, implying 
10 km, 1.4 Mq NS; this tests for consistency with the standard model. 
a These parameters vary among W37's observations, and are presented in Table l3l 
b This parameter reached the hard limit of the model. 
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TABLE 3 
Additional Parameters for W37 



Obs. 


Cts 


F x , (0.5-2.5 keV) 


N H 


Intervals of spectral extractions 






(1CT 14 orgs s- 1 ) 


(10 20 cm" 2 ) 


(10 3 s) 



29.006:43.006, 53.006:66.102 
0:29.006, 43.006:53.006 
74.104:140.209 
266.385:278.006, 287.506:289.506 
278.006:287.506, 289.506:332.490 
1038.010:1045.010 
980.879:1038.010, 1045.010:1051.550 



Note. — Spectral fits to various parts of W37 data, using Llovd (2003) hydrogen-atmosphere 
neutron star model and power-law model (as in Table 121 . with Nfj the only variable parameter. 
Times are barycentered and should be added to reference time 149706994.3 (seconds after MJD 
50814.0). All errors are 90% confidence limits. Distance of 4.85 kpc is assumed. 
a This parameter reached the hard limit of the model. 
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